MATHEMATICAL ENGINEERING TECHNICAL REPORTS Estimation of Exponential-Polynomial Distribution by Holonomic Gradient Descent
نویسندگان
چکیده
We study holonomic gradient decent for maximum likelihood estimation of exponentialpolynomial distribution, whose density is the exponential function of a polynomial in the random variable. We first consider the case that the support of the distribution is the set of positive reals. We show that the maximum likelihood estimate (MLE) can be easily computed by the holonomic gradient descent, even though the normalizing constant of this family does not have a closed-form expression and discuss determination of the degree of the polynomial based on the score test statistic. Then we present extensions to the whole real line and to the bivariate distribution on the positive orthant.
منابع مشابه
MATHEMATICAL ENGINEERING TECHNICAL REPORTS Holonomic Gradient Method for the Distribution Function of the Largest Root of a Wishart Matrix
We apply the holonomic gradient method introduced by Nakayama et al. [19] to the evaluation of the exact distribution function of the largest root of a Wishart matrix, which involves the hypergeometric function 1F1 of a matrix argument. Compared to methods based on infinite series expansion in terms of zonal polynomials, whose convergence is slow, the proposed method is more practical for compu...
متن کاملMATHEMATICAL ENGINEERING TECHNICAL REPORTS Holonomic Gradient Method for Distribution Function of a Weighted Sum of Noncentral Chi-square Random Variables
We apply the holonomic gradient method to compute the distribution function of a weighted sum of independent noncentral chi-square random variables. It is the distribution function of the squared length of a multivariate normal random vector. We treat this distribution as an integral of the normalizing constant of the Fisher-Bingham distribution on the unit sphere and make use of the partial di...
متن کاملMATHEMATICAL ENGINEERING TECHNICAL REPORTS Calculation of Orthant Probabilities by the Holonomic Gradient Method
We apply the holonomic gradient method (HGM) introduced by [9] to the calculation of orthant probabilities of multivariate normal distribution. The holonomic gradient method applied to orthant probabilities is found to be a variant of Plackett’s recurrence relation ([14]). However a straightforward implementation of the method yields recurrence relations more suitable for numerical computation ...
متن کاملAn eigenvalue study on the sufficient descent property of a modified Polak-Ribière-Polyak conjugate gradient method
Based on an eigenvalue analysis, a new proof for the sufficient descent property of the modified Polak-Ribière-Polyak conjugate gradient method proposed by Yu et al. is presented.
متن کاملImproving the Performance of Bayesian Estimation Methods in Estimations of Shift Point and Comparison with MLE Approach
A Bayesian analysis is used to detect a change-point in a sequence of independent random variables from exponential distributions. In This paper, we try to estimate change point which occurs in any sequence of independent exponential observations. The Bayes estimators are derived for change point, the rate of exponential distribution before shift and the rate of exponential distribution after s...
متن کامل